Differential Effects of Senescence on the Phloem Exports of Cadmium and Zinc from Leaves to Grains in Rice during Grain Filling

In rice, non-essential toxic cadmium (Cd) and the essential nutrient zinc (Zn) share similar transport pathways, which makes it challenging to differentially regulate the allocation of these elements to the grain. The phloem is the main pathway for the loading of these elements into rice grains. It has long been accepted that tissue senescence makes the nutrients (e.g., Zn) stored in leaves available for further phloem export toward the grain. Whether senescence could drive the phloem export of Cd remains unclear. To this end, the stable isotopes 111Cd and 67Zn were used to trace the phloem export and the subsequent allocation of Cd and Zn from the flag leaves, where senescence was accelerated by spraying abscisic acid. Furthermore, changes upon senescence in the distribution of these elements among the leaf subcellular fractions and in the expression of key transporter genes were investigated. Abscisic acid-induced senescence enhanced the phloem export of Zn but had no impact on that of Cd, which was explained by the significant release of Zn from the chloroplast and cytosol fractions (concentrations decreased by ~50%) but a strong allocation of Cd to the cell wall fraction (concentration increased by ~90%) during senescence. Nevertheless, neither Zn nor Cd concentrations in the grain were affected, since senescence strengthened the sequestration of phloem-exported Zn in the uppermost node, but did not impact that of phloem-exported Cd. This study suggests that the agronomic strategies affecting tissue senescence could be utilized to differentially regulate Cd and Zn allocation in rice during grain filling.


Introduction
Mineral levels in cereals are critical for food safety and nutrition, among which the levels of cadmium (Cd) and zinc (Zn) are of particular concern [1,2]. For both higher plants and humans, Cd is a non-essential and toxic metal that could cause metabolic disorders, whereas Zn is an essential nutrient that is fundamental to enzyme activities [3,4]. Among the world's population, especially the Asian population, rice accounts for a major portion of the caloric intake [5]. However, rice is characterized as a major source of human exposure to Cd but a poor source of dietary Zn [2,6]. Hence, Cd accumulation must be diminished and Zn accumulation be fortified in rice grains, in order to prevent the population from experiencing Cd toxicities and Zn deficiencies. This requires a further understanding of the allocation processes of these elements in rice plants.
Owing to their similar chemical properties, Cd and Zn transport processes from soil to grains are generally similar, resulting in a simultaneous reduction of Zn in grains when agronomic measures are taken to reduce Cd in grains [7][8][9]. Thus far, it remains challenging to differentially regulate Cd and Zn accumulations in rice grains [10]. In acidic soils, both the solubility of Cd and Zn and the subsequent uptake decrease with an increase in soil pH, since the higher pH causes weaker ion-proton exchange and thereby stronger and nonspecific adsorption of Cd and Zn onto the negatively charged solid phase [7,11]. In plants, it is well known that Cd shares transporters with Zn, iron (Fe) and manganese [6]. For example, it was found that both Zn and Cd can be taken up by OsIRT1 in rice roots [12], where OsHMA3 mediates vacuole compartmentalization while OsHMA2 is involved in the xylem loading of both elements for root-to-shoot translocation [1,13]. Both Cd and Zn are translocated in the xylem mostly as free ions [14,15] and highly allocated toward the photosynthetic organs (e.g., leaves) with the transpiration flow during the vegetative growth period [16,17]. During grain filling, OsHMA2 in the node facilitates the allocation of both Cd and Zn to the grain by mediating the intervascular transport of these elements [18]. Together, these similarities make the Cd and Zn content in grains always positively correlated. Indeed, some studies found that applying lime reduces the soilavailable Cd and grain Cd concentrations, but also reduces those of Zn [7,9]. To produce grains with low Cd and high Zn levels, it is therefore necessary to reveal the physiological mechanisms that could differentially regulate the allocation of Cd and Zn to grains.
Relatively less attention has been paid to the phloem transport of Cd and Zn, despite the phloem being the main pathway for the loading of these elements into rice grains [19,20]. Previous studies have shown that the majority (>90%) of Cd and Zn in rice grains is allocated via the phloem after having been exported from phloem source organs or after xylem-to-phloem transfer in the nodes [21,22]. In rice during grain filling, the leaves are the primary phloem source that significantly exports Cd and Zn to the phloem sink, i.e., grains [23,24]. To date, it remains unclear if there are any differences between Cd and Zn in their phloem transport that could be used for differentially regulating the allocation of these elements during grain filling.
It is well accepted that tissue senescence makes the minerals stored in leaves available for further phloem export [25]. During leaf senescence, autophagy catabolizes cellular components gradually [26]. In general, the chloroplast is the first target of autophagic degradation, which turns the leaf from green to yellow [25]. Thus, the minerals intensively contained in the chloroplasts are the first to become available for phloem transport. For example, more than half of nitrogen in leaves is present in Rubisco in chloroplasts, in which phloem export is well synchronized with the onset of leaf senescence [27]. With respect to Zn, the dependency between its phloem export and tissue senescence has long been recognized in cereal crops [28,29]. Indeed, accelerating leaf senescence is suggested as an effective strategy to improve the grain nutrition value for protein, Fe and Zn [30]. By contrast, Cd is relatively less mobile for phloem export from leaves. In rice, Wiggenhauser, et al. [31] revealed that Cd can complex with organic acids in the vacuole or with O-donors in the cell wall, thereby contributing to the immobility of Cd in leaves and stems. These cellular components are degraded or even remain intact at the late stage of autophagy [26], which may lead to the immobility of Cd during senescence. Upon the phloem export from the leaves, the nodes, particularly the uppermost node I, allocate Cd and Zn further to the grain or other organs. Although this allocation process is mainly regulated by OsHMA2, which is shared by Cd and Zn as mentioned above, there are transporters that are highly specific for Cd (OsLCT1) and for Zn (OsZIP3) in regulating the intervascular transfer of these elements [32,33]. When the phloem export of minerals changes during tissue senescence, plants are able to coordinate the expression of nodal transporter genes to meet the sink demand or avoid potential toxicity [34]. Taken together, we hypothesized that the difference between the phloem transport of Cd and Zn relies on their distinct responses to senescence during grain filling.
Plants 2023, 12, 1902 3 of 16 This study aimed at (1) testing whether Cd and Zn respond differently to senescence in terms of phloem export and subsequent allocation, and (2) understanding the physiological processes that lead to the differences. To this end, the stable isotopes 111 Cd and 67 Zn were used to trace the allocation of Cd and Zn from the flag leaf to other parts of rice. Plant senescence was accelerated by applying the phytohormone abscisic acid. The subcellular distribution of Cd and Zn tracers in the flag leaf and the expression of some key transporter genes were also investigated upon tissue senescence.

Plant Growth and Total Cd and Zn Accumulation
During the entire growth period, the plant grew healthily and no poisoning symptoms of Cd were observed. From the 13th day after anthesis, the relative chlorophyll concentration (SPAD unit) of labeled flag leaves was significantly different between the ABA treatment and the blank control (CK; Figure 1). The SPAD unit decreased by 48% in the ABA plants during the 9 days of hormone treatment, but by only 19% in the CK. that the difference between the phloem transport of Cd and Zn relies on their distinct responses to senescence during grain filling.
This study aimed at (1) testing whether Cd and Zn respond differently to senescence in terms of phloem export and subsequent allocation, and (2) understanding the physiological processes that lead to the differences. To this end, the stable isotopes 111 Cd and 67 Zn were used to trace the allocation of Cd and Zn from the flag leaf to other parts of rice. Plant senescence was accelerated by applying the phytohormone abscisic acid. The subcellular distribution of Cd and Zn tracers in the flag leaf and the expression of some key transporter genes were also investigated upon tissue senescence.

Plant Growth and Total Cd and Zn Accumulation
During the entire growth period, the plant grew healthily and no poisoning symptoms of Cd were observed. From the 13th day after anthesis, the relative chlorophyll concentration (SPAD unit) of labeled flag leaves was significantly different between the ABA treatment and the blank control (CK; Figure 1). The SPAD unit decreased by 48% in the ABA plants during the 9 days of hormone treatment, but by only 19% in the CK. There were no significant differences in the grain yield and plant biomass between ABA and CK treatments ( Table 1). The ABA treatment had no impact on the distribution of dry matter among plant organs. Biomass was mostly allocated to the grains in the labeling tiller (including grains, H + R, flag leaves, node-I and stems). For the whole plant, over 75% of biomass was distributed in the shoots (8.81-8.91 g). There were no significant differences in the grain yield and plant biomass between ABA and CK treatments ( Table 1). The ABA treatment had no impact on the distribution of dry matter among plant organs. Biomass was mostly allocated to the grains in the labeling tiller (including grains, H + R, flag leaves, node-I and stems). For the whole plant, over 75% of biomass was distributed in the shoots (8.81-8.91 g). In the two treatments, except for the flag leaves labeled with Cd, the total Cd concentration was the highest in the root ( Table 2). The total Zn concentration was the highest in node-I, followed by the root. Between 6 and 15 DAA, the total amount of Cd and Zn in grains kept increasing, while no significant differences were observed between the treatments (Table S1). By contrast, compared to the CK, the total amounts of Zn were increased in the node-I by 44% and in the stem by 63% upon the ABA treatment, whereas the amounts of Cd in these organs were not affected. In addition, while the concentrations of Cd were not affected by the ABA, the Zn concentrations were increased by 56% in the other tillers. In addition, in the other tillers and node-I, while the Cd concentrations were not affected by the ABA treatment, the Zn concentrations significantly increased by 56% and 39%, respectively. After 24 h of isotope-labeled Cd (Cd lab ) and labeled Zn (Zn lab ) application, 0.82 µg Cd lab and 1.20 µg Zn lab were absorbed (Table S2).

Phloem Export and Allocation of Foliar Applied Cd and Zn
In this experiment, the concentration of Cd lab and Zn lab in the phloem exudates was standardized by Fe. The Cd lab concentration in the phloem exudate was similar between the treatments ( Figure 2). By contrast, compared to CK, the Zn lab concentration in the phloem exudate increased by 110% under ABA treatment (p < 0.05). After 24 h of isotope-labeled Cd (Cdlab) and labeled Zn (Znlab) application, 0.82 µg Cdlab and 1.20 µg Znlab were absorbed (Table S2).

Phloem Export and Allocation of Foliar Applied Cd and Zn
In this experiment, the concentration of Cdlab and Znlab in the phloem exudates was standardized by Fe. The Cdlab concentration in the phloem exudate was similar between the treatments ( Figure 2). By contrast, compared to CK, the Znlab concentration in the phloem exudate increased by 110% under ABA treatment (p < 0.05).  Tables 3 and S3, there were no significant differences between the treatments in the amount of Cdlab (0.10-0.12 µg) and Znlab (0.81-0.82 µg) remobilized from the flag leaves. The Znlab exported from flag leaves was distributed into various parts of the rice. The Znlab concentrations were 0.6-0.8 µg kg −1 in the rhizosphere soil where, however, no Cdlab was detected. At maturity, the exported Znlab was mainly found in the grains, stems, node-I, other tillers and roots, whereas Cdlab was found only in the grains, H + R and stems ( Figure 3). In the ABA treatment, the amount of Znlab in node-I was significantly higher than that that in the CK (Table 3). In the ABA treatment, 57% of the exported Znlab was allocated to the node-I was significantly higher than the fraction (29%) in the CK (Figure 3b).  Tables 3 and S3, there were no significant differences between the treatments in the amount of Cd lab (0.10-0.12 µg) and Zn lab (0.81-0.82 µg) remobilized from the flag leaves. The Zn lab exported from flag leaves was distributed into various parts of the rice. The Zn lab concentrations were 0.6-0.8 µg kg −1 in the rhizosphere soil where, however, no Cd lab was detected. At maturity, the exported Zn lab was mainly found in the grains, stems, node-I, other tillers and roots, whereas Cd lab was found only in the grains, H + R and stems ( Figure 3). In the ABA treatment, the amount of Zn lab in node-I was significantly higher than that that in the CK (Table 3). In the ABA treatment, 57% of the exported Zn lab was allocated to the node-I was significantly higher than the fraction (29%) in the CK (Figure 3b).

Subcellular Distribution of Foliar-Applied Cd and Zn
The ABA treatment significantly affected the concentrations of Cdlab and Znlab in the subcellular fractions extracted from the flag leaves ( Figure 4). In the cell wall fraction, the concentration of Cdlab increased by 89% but that of Znlab remained stable with ABA application compared to CK. Regardless of the treatment, more than 50% of Cdlab was sequestrated in the cell wall ( Figure S1). By contrast, in the chloroplast and cytosol fractions, ABA treatment significantly lowered the concentration of Znlab by more than 50% but had no

Subcellular Distribution of Foliar-Applied Cd and Zn
The ABA treatment significantly affected the concentrations of Cd lab and Zn lab in the subcellular fractions extracted from the flag leaves ( Figure 4). In the cell wall fraction, the concentration of Cd lab increased by 89% but that of Zn lab remained stable with ABA application compared to CK. Regardless of the treatment, more than 50% of Cd lab was sequestrated in the cell wall ( Figure S1). By contrast, in the chloroplast and cytosol fractions, ABA treatment significantly lowered the concentration of Zn lab by more than 50% but had no impact on the level of Cd lab . In the organelle fraction, for both Cd lab and Zn lab , the concentrations were similar between the treatments. impact on the level of Cdlab. In the organelle fraction, for both Cdlab and Znlab, the con trations were similar between the treatments.

Expression Profile of OsHMA2, OsLCT1 and OsZIP3 in the Flag Leaves and Node-I
In the flag leaves, the ABA treatment significantly increased the relative expres level of OsHMA2 by 2.2 times compared to the CK, but had no impact on the expres of OsLCT1 (Figure 5a,b). In the node-I, both the expression levels of OsHMA2 and OsL were promoted under the ABA treatment by around two times (Figure 5c,d). The rela expression level of OsZIP3 remained stable regardless of the treatments (Figure 5e).

Expression Profile of OsHMA2, OsLCT1 and OsZIP3 in the Flag Leaves and Node-I
In the flag leaves, the ABA treatment significantly increased the relative expression level of OsHMA2 by 2.2 times compared to the CK, but had no impact on the expression of OsLCT1 (Figure 5a,b). In the node-I, both the expression levels of OsHMA2 and OsLCT1 were promoted under the ABA treatment by around two times (Figure 5c,d). The relative expression level of OsZIP3 remained stable regardless of the treatments (Figure 5e).

Plant Culture
Rice (Oryza sativa L.) cultivar 'Meixiangzhan 2′ was used in this study, which has a higher risk of accumulating Cd in the grain to levels exceeding the regulatory limit (Li

Plant Culture
Rice (Oryza sativa L.) cultivar 'Meixiangzhan 2' was used in this study, which has a higher risk of accumulating Cd in the grain to levels exceeding the regulatory limit (Li Guijie, 2016). Seeds were surface sterilized by using 2% NaClO solution for 20 min and rinsed five times with deionized water. The seeds were germinated in the dark for The seedlings were transplanted at the tillering stage (33 days after germination) to pots (155 mm in height and 162 mm in diameter; one plant per pot) filled with substrate mixing three volumes of peat substrate (Fine, Jiffy, Zwijndrecht, Holland) and one volume of medium grain quartz sand (particle size 0.5-1 mm). Slow-release fertilizer Plantacote Depote 6M (NPK 14-9-15, Meyer, Rellingen, Germany) was supplied to the substrate at the rate of 1 g per liter, and 0.1 g per liter micronutrient fertilizer Fetrilon-Combi (0.7% Mg, 1.5% B, 0.5% Cu, 4% Fe, 3% Mn, 0.05% Mo, 4% Zn, COMPO EXPERT GmbH, Poznan, Poland) was supplied. The rice plants were planted in an outdoor plastic greenhouse from April to early July 2022. The rice was exposed to natural sunlight and ventilation in Guangzhou (23 • 26 N, 113 • 23 E). The daily average temperature was 26.8-31.0 • C (min 15 • C, max 38 • C) during cultivation. The monthly average relative humidity was 75%-80% (data: Chinese National Meteorological Centre). The rice was re-flooded every dusk with tap water.
The phenology was carefully observed for each plant. We defined 'anthesis' when the anthers were visible, and the tiller carrying the first anthesis floret was defined as the main tiller [36]. When 3 florets carried anthers, these 3 tillers were kept while other tillers were removed by cutting.

Stable Cadmium and Zinc Isotope Labeling
The 111 Cd-enriched isotope and 67 Zn-enriched isotope were purchased from Isoflex (the United States). The isotopic abundances of Cd and Zn are shown in Table S4. A solution containing 1.0 mM 111 Cd and 2.0 mM 67 Zn was prepared and adjusted to neutral (pH = 6.5) using MES and KOH. To enhance the surface adhesion of solutions toward leaves, Tween 80 (0.05% v/v) was added. At 5 days after anthesis (DAA), ten 5 µL droplets of the isotope solution were applied to the flag leaves of main tillers. Note that the other 2 tillers were not labeled. This time was selected considering that grain accumulates dry matter rapidly between 4 and 18 DAA [37,38]. After 24 h, the treated flag leaves were rinsed with a solution containing 2% nitric acid (HNO 3 ) and 3% ethanol to remove the unabsorbed isotope label [39,40].

Plant Hormone Treatment
Two treatments were established with 10 replicates each. Abscisic acid (ABA) was used to accelerate the plant senescence [41]. ABA solution was prepared by dissolving dextrorotatory ABA powder in NaOH and diluted to 200 µM with ultrapure water. Meanwhile, a solution with the same concentration of NaOH was prepared as the blank control (CK). These solutions also contained Tween 80 (0.05% v/v), which can reduce the surface tension of rice leaves effectively, and prolong the residence time of hormones on the leaf surface. For all replicates of each treatment, in total 250 mL of the solution was sprayed on the plants from the top side every 24 h, which lasted 9 days. Before the hormone treatment, the chlorophyll relative concentration of each plant was measured using a soil and plant analyzer development (SPAD-502, Minolta, Milton Keynes, UK). The flag leaves of the plants were taken as the measurement objects, and 20 points were randomly measured and averaged, each group was measured five times and recorded. The hormone treatment was continued for 9 days, until the SPAD values were sufficiently distinguished between the treatments.

Phloem Exudate Extraction
Flag leaves were collected freshly after the hormone treatment (at 15 DAA) for phloem exudate extraction, according to Tetyuk, et al. [42]. Briefly, for each treatment, 4 pieces of fresh flag leaves (including the blade and sheath) from the main tillers of different plants subjected to the same treatment were cut at the node level and recut (about 1 mm) at the sheath bottom in a petri dish containing 20 mM K 2 -EDTA solution. Then, the leaves were quickly transferred to a centrifuge tube (5 mL, Eppendorf, Hamburg, Germany) containing 1.4 mL of 20 mM K 2 -EDTA. After immersing the sheath bottom in the solution, the leaves were placed in a dark and humid environment for 1 h. After bath washing the sheath bottom using flowing ultrapure water, the flag leaves were transferred to a new tube containing 1.4 mL ultrapure water. The phloem exudate was extracted to the water for 5 h in a dark and humid environment. Deng, et al. [43] suggested that Fe can be used as a reference element to standardize the variations in the extraction efficiency as the concentration of Fe is relative constant in the phloem exudates. Therefore, in the present study, the concentrations of Cd lab and Zn lab in the phloem exudates were standardized and expressed as the ratios between Zn/Cd and Fe concentrations.

Extraction of Subcellular Fractions
From the extras' replicates, the flag leaves were collected freshly for the extraction of subcellular fractions, according to the method used in Wu, et al. [44]. Briefly, each fresh flag leaf was homogenized with 10.0 mL of precooled (4 • C) extraction solution containing 250 mmol·L −1 sucrose, 50 mmol·L −1 Tris HCl (pH = 7.50) and 1 mmol·L −1 dithiocarbaryl. The homogenate was filtered by a 240 µM nylon cloth, the rest in the nylon cloth was the cell wall fraction that may contain the cell walls and cell wall debris. After that, the filtrate was collected into a 50 mL centrifuge tube, centrifuged at 1500× g with a high-speed freezing centrifuge for 10 min, and the precipitation was the chloroplast fraction. The supernatant was further centrifuged at 15,000× g for 30 min. The precipitation was the organelle fraction that may contain the cell membrane and organelles (apart from the chloroplast). The supernatant was the cytosol fraction that may contain the cytoplasm and vacuoles. All steps were performed at 4 • C.

Plant Harvesting
To reveal the allocation of Cd and Zn upon the hormone treatments, the plants were harvested at 6 DAA (after anthesis) and at 15 DAA (after anthesis). The roots were carefully separated from the soil with a tap-water flow and deionized water, then bathed in 4 L of 1 µM DTPA solution for 10 min and rinsed with deionized water for 1-2 min to remove the Cd adsorbed onto the apoplast [45]. For the shoots, the labeling tillers were divided into grains, husk + rachis (H + R), flag leaves, node-I (the uppermost node), stems (including the lower nodes and internodes), and the lower leaves, while the other tillers were mixed as one sample. All these shoot organs were washed with tap water followed by deionized water. Then, all the plant samples were dried to constant weight in an oven at 70 • C, weighed, and ground with tungsten carbide cups and balls. In addition, rhizosphere soil was collected within 2 cm surrounding the rice root [46], air-dried, and sieved to <80 µm.

Sample Processing
The samples of phloem exudate were diluted to 5 mL with deionized water, and 100 µL concentrated HNO 3 was added. The phloem exudate samples were determined for the total concentration of Zn, Fe and Cd, and for the isotope ratio of 111 Cd/ 110 Cd and 67 Zn/ 66 Zn.
The subcellular fractions samples were digested with 5 mL nitric acid by an electric heating plate (JinRongYuan, JRY-D350-D, Changsha, China) at 105 • C for 180 min and dilute to 20 mL to determine the total concentration of Zn and Cd, and the isotope ratio of 111 Cd/ 110 Cd and 67 Zn/ 66 Zn. About 0.10 g of the plant samples and soil samples were preliminarily digested for >12 h in 1 mL 69% HNO 3 and 3 mL ultrapure water. After preliminary digestion, 1 mL 30% (w/v) H 2 O 2 was added. The samples were digested by a microwave digestion extraction system (CEM, MARS6, Boston, MA, USA) at 120 • C for 10 min, 180 • C for 10 min and 210 • C for 20 min. The digests were diluted to 20 mL using ultrapure water and filtrated through water phase needle filters (0.45 µm).

RNA Extraction and qRT-PCR Analysis
The RNA extraction and qRT-PCR analysis were carried out according to Sun, et al. [47]. Total RNA was extracted from fresh plant node-I and flag leaves using the FastPure Plant Total RNA Isolation Kit (Vazyme, Nanjing, China) and then reverse-transcribed using the HiScript III RT SuperMix for qPCR (Vazyme, Nanjing, China). qRT-PCR was performed using the StepOne Real-Time PCR System (StepOne, Shanghai, China). The relative expression levels were computed by the 2 −∆∆CT method of relative quantification. Oryza sativa L. actin gene (OsActin) was used as the internal control, with primers being shown in Table S5.

Data Analysis
For each organ of each plant, the amount of total Cd (Q where A refers to the abundance of the isotope provided in Table S4.  (6) For Cd lab and Zn lab the net amount exported (Q exp Cd and Q exp Zn ) from the labeled flag leaf was estimated as the decreased amount of Cd lab /Zn lab in the flag leaf before and after the hormone treatment (i.e., between 6 and 15 DAA), given that the decrease was statistically significant.
Statistical analysis with Student's t-test was used to identify the significant differences between the treatments. Statistical analyses were conducted using Excel 2016 (Microsoft, Washington, DC, USA) for Windows, and all graphs were drawn by OriginPro 2023 (Origin-Lab, Northampton, MA, USA).

ABA-Induced Senescence Fortified the Zn Phloem Export, but Had No Impact on the Cd Phloem Export
Compared to CK, the significantly decreased SPAD value showed that ABA significantly accelerated the flag leaves' senescence ( Figure 1). Previous studies showed that senescence promotes the remobilization of minerals from leaves in crops [33,48,49]. Phloem is the only pathway of Zn and Cd remobilized from leaves [22,50]. In the present study, the ABA treatment markedly enhanced the phloem export of Zn lab applied to the flag leaves, but had no impact on the phloem export of Cd lab (Figure 2). The changes in the subcellular distribution of Cd and Zn upon ABA-induced senescence were investigated to understand the physiological mechanism behind this observation.
It is generally proposed that Zn is mainly present in the cytosol and functional organelles (e.g., chloroplast and mitochondrion), given its essentiality for enzyme functioning, while the potentially toxic Cd is mainly immobilized by the cell wall or compartmentalized in the vacuole [3,51]. ABA-induced senescence promoted the release of Zn lab in chloroplasts and cytosol components significantly. The total proportion of Zn lab allocated to these components in the CK was 64%, while ABA-induced senescence resulted in a large outflow of Zn lab from both of them ( Figure 4) and a significant decrease in the proportion (down to 22%; Figure S1). Consequently, the phloem export of Zn lab markedly increased ( Figure 2). By contrast, ABA-induced senescence promoted the distribution and fixation of Cd lab in flag leaves to the cell wall component, which may adsorb some, if not all, Cd lab released by the other cellular components during senescence, thus limiting the phloem export of Cd lab . The ABA treatment promoted the release of Cd from some subcellular components, as the concentrations of Cd lab in the chloroplast and cytosol components decreased (but the significance may be diluted by variability among replicates) (Figure 4b,d). In the flag leaves, Cd lab was mainly stored in the cell wall component (Figure 4a), which was similar to many other studies [52,53]. In addition, cell wall precipitation is a defense mechanism of plant cells against heavy metal toxicity [54]. During grain filling, with the senescence of rice, the Cd storage as Cd-O in the cell wall may contribute to the immobilization of Cd in the flag leaves [31]. The remobilization of Zn in flag leaves may depend on the plant's developmental stage [55]. Similar transportation of Zn and Cd in the phloem has been shown by Page and Feller [56]. It has been found that Cd in the older leaves was subsequently remobilized in the phloem only at trace levels compared to Zn [56]. This study suggests that ABA-induced senescence did not promote the release of Cd from leaf cells, which may explain the limited phloem export of Cd. Therefore, compared to Cd, Zn appeared to be more mobile in the flag leaves during ABA-induced senescence, which may explain the differential effect of senescence on the phloem export of Cd and Zn from the flag leaves.
However, the enhanced Zn lab phloem export did not result in more Zn lab export from the flag leaves (Table 3). This may be due to a balanced import and export of Zn lab mass in the flag leaves. In this study, we found that Zn lab applied to the flag leaf was allocated to the roots, which was similar to Haslett [57] reported in wheat. Meanwhile, according to Table S1, compared with CK, the root-to-shoot ratio of total Zn amount increased significantly in the ABA treatment, indicating that more Zn lab was translocated from the root, which may result in the absence of apparent Zn lab export in the flag leaves. Therefore, the mass balance data cannot distinguish between the leaf source and root source, but the Zn lab in the phloem sap can be determined to be exported from the leaves.

ABA-Induced Senescence Affects the Allocation of Phloem-Exported Zn by Altering the Expression of Transporter Genes in Flag Leaves and Nodes
In accordance with previous studies, Cd and Zn were found to be highly enriched in the uppermost Node-I compared to other organs [21,58,59]. The node-I is known as a key hub for the allocation of Cd and Zn in rice. In addition, the flag leaves were the phloem sources of Cd and Zn in the present study. Given these, in the node-I and flag leaves, we investigated the effect of ABA-induced senescence on the relative expression levels of some key transporter genes (OsLCT1, OsHMA2 and OsZIP3).
OsLCT1 is known to be mainly expressed in leaf blades and nodes and facilitates the export of Cd, but not Zn [13,32]. The expression level for OsLCT1 in the flag leaves was not affected by ABA-induced senescence (Figure 5b), which corresponded to the observation that senescence had no effect on Cd export from the flag leaves found in our study ( Figure 2). ABA-induced senescence promoted the expression of OsLCT1 in the node-I (Figure 5d), but had no effect on the distribution of Cd in node-I (Table 3). It was reported that decreased expression of OsLCT1 can lead to a decrease in the concentration of Cd in the phloem, and OsLCT1 in the diffuse vascular bundles of node-I is important for Cd transport into grains [32]. We found that senescence did not promote the expression of OsLCT1 in flag leaves (Figure 5b). This result may explain why ABA-induced senescence had no effect on the phloem export of Cd lab from the flag leaves and the allocation of Cd lab through the phloem toward the other organs (e.g., node-I, grains).
Although ABA-induced senescence promoted the export of Zn in the phloem, we did not find a significant increase in Zn in the grains. This is mainly due to the significant accumulation of Zn in the node-I in the ABA treatment (Table 3). OsHMA2 in flag leaves was strongly upregulated by ABA-induced senescence (Figure 5a). OsHMA2 has been proven to play a key role in the phloem transport of Zn and Cd, expressed mainly in node-I at the flowering stage, and also expressed in flag leaves [13]. OsHMA2 is localized at the enlarged and diffuse vascular bundles in the nodes and regulates the preferential distribution of Cd and Zn through the phloem to the developing tissues [18]. It was found that the root-to-shoot translocation of Zn via OsHMA2 remained stable when the concentration of Cd in the nutrient solution increased from 0 to 4 nM, indicating that Cd can hardly compete with Zn for OsHMA2 at far lower concentrations than Zn [60]. In this study, the concentration of total Zn in flag leaves was much higher than that of total Cd (The ratio of total Zn to total Cd concentration was greater than 10). Therefore, compared to Cd, Zn was more competitive for the binding site of OsHMA2, which may mainly promote the transport of Zn in the phloem.
ABA-induced senescence had no effect on the expression level of OsZIP3 in the node-I (Figure 5e). Sasaki, et al. [33] found that OsZIP3 is localized in rice nodes and is an essential transporter responsible for the distribution of Zn to developing tissues through unloading Zn from the xylem. In this study, while ABA treatment significantly promoted the rootto-shoot translocation of Zn, it did not increase the level of zinc in grains (Table S1). This might be due to the unaffected expression level of ZIP3 in node-I, limiting the unloading Zn from the xylem the subsequent allocation of Zn toward grains.

Conclusions
Senescence induced by abscisic acid (ABA) promoted the phloem export of Zn from flag leaves, but not Cd, because senescence caused significant Zn release from the chloroplast and cytosol fractions of leaf cells, while Cd was immobilized in the cell wall fractions. Meanwhile, the expression of OsHMA2 in the flag leaves was promoted by senescence, which facilitated the loading of Zn into the phloem. Nevertheless, due to the strong se-questration of phloem-exported Zn in the uppermost node-I, the accumulation of Zn in the grain was not affected by accelerated senescence. By contrast, although senescence enhanced OsLCT1 expression in node-I, the allocation of phloem-exported Cd was also not different in the node-I thanks to the limiting OsLCT1 expression of flag leaves upon ABA treatment. These results showed that agronomic measures, such as irrigation and fertilization, could be used to differentially regulate the allocation of Zn and Cd in rice during grain filling through their effects on plant senescence. To effectively manage the grain Zn and Cd levels, however, further studies are still required to better understand the response of node-I in distributing the phloem-exported Zn and Cd toward plant organs.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/plants12091902/s1, Table S1: Amount of total Cd and Zn in each organ harvested at 6 DAA (after anthesis) and at 15 DAA (after anthesis) in rice; Table S2: Cd lab and Zn lab absorption in flag leaves during 24 h; Table S3: The amount Cd lab and Zn lab among the rice in two treatments at 6 DAA (after anthesis) and at 15 DAA (after anthesis); Table S4: isotopic abundances of Zn and Cd (data from the isotope certificate); Table S5: Primers for qRT-PCR. Figure S1: Subcellular distribution proportion of Cd lab (a) and Zn lab (b) in flag leaves of rice at 15DAA (after anthesis).

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.